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Abstract 

We study characteristic features of the eigenvalues of the Wilson- 
Dirac operator in topologically non-trivial gauge field configurations 
by examining complete spectra of the fermion matrix. In particular 
we discuss the role of eigenvectors with real eigenvalues as the lattice 
equivalents of the continuum zero-modes. We demonstrate, that those 
properties of the spectrum which correspond to non-trivial topology 
are stable under adding fluctuations to the gauge fields. The behavior 
of the spectrum in a fully quantized theory is discussed using QED2 
as an example. 
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1 Introduction 



The Atiyah-Singer index theorem |Q is considered to be among the greatest 
achievements of differential geometry. It has entered the physics hterature 
in many places. Of particular interest is its application in gauge theories 
with fermions. The Atiyah-Singer index theorem relates the topological 
charge I'lA] of a sufficiently smooth gauge field configuration A to the index 
of the Dirac operator on some compact manifold (or for gauge fields with 
suitable boundary conditions). The index of the operator can be expressed 
as the difference of the number of (lineraly independent) negative (n_) and 
positive (n+) chirality zero modes. The zero modes are eigenstates of the 
Dirac operator with eigenvalue E = 0. Since 75 anti-commutes with ij^, the 
zero modes can be chosen as eigenstates V'+j ^- of 75 with 75^^+ = +tp+ and 
75 ^/^_ = n+ (n_) denotes the number of independent states -0+ (V'-)- 

The Atiyah-Singer index theorem Q then reads 

u[A] = n_ - n+ . (1) 

This relation allows to relate properties of fermionic observables to the topol- 
ogy of the background gauge field. Unfortunately, the gauge fields carry- 
ing the measure of the continuum path integral do not obey the necessary 
smoothnes condition for the index theorem to apply in the fully quan- 
tized theory. The application of the index theorem is thus restricted to 
semiclassical arguments. 

On the lattice the situation is different. It is possible to assign a topolog- 
ical charge to all gauge field configurations (see e.g. for a recent review), 
up to so-called exceptional configurations which do not contribute in the 
continuum limit. The lattice path integral for the gauge fields thus effec- 
tively can be decomposed into topological sectors. Certainly it would be 
valuable to be able to understand the behavior of fermionic observables in 
each sector. Unfortunately there is no analytic result for an index theorem 
on the lattice. However, many contributions can be found in the literature 
IQ] - ||2^, where a probabilistic manifestation of the Atiyah-Singer index 
theorem on the lattice and its consequences are discussed. 

The pioneering papers ^, ^ investigated the manifestation of the index 
theorem on the lattice by analyzing the pseudoscalar density operator which 
through chiral Ward identities can be related to the topological charge. In 
two interesting papers by Smit and Vink these relations were extended 
to a lattice version of the Witten-Veneziano formula for the pseudoscalar 
masses. In two numerical studies the interplay between pseudoscalar 
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density operators and the topological charge was further investigated for 
two-dimensional models. Complete spectra of the Wilson-Dirac operator 
for non-abelian gauge fields first appeared in 

An important step in understanding the manifestation of the index the- 
orem on the lattice was to realize ^ |lO| that the eigenvectors of the 
lattice Dirac operator with real eigenvalues should be interpreted as the 
lattice counterparts of the continuum zero modes. This result can best be 
understood [l^ by using the fact that for the Wilson-Dirac operator only 
eigenvectors -0 with real eigenvalues can have non-vanishing pseudoscalar 
matrix elements (■0,75V')! similar to the zero modes in the continuum. It 
was furthermore realized ||5|, |l^, that only real eigenvalues of D can give 
rise to exact zero eigenvalues for properly chosen bare mass. As is explicitly 



shown in [|10[ this can be used to compute the topological charge by counting 
the level crossings of the eigenvalues of the numerically simpler hermitian 
operator 75Z) as the bare quark mass is varied. This method was later re- 



assessed in the framework of the overlap formalism |11|. The procedure was 
successfully applied to study the interplay between topological charge and 
the spectrum of 75D in SU(A^) gauge theories in various settings including 
2d models ]12| as well as static and thermalized configurations for SU(2) and 
SU(3) In the last reference in |11| the index theorem is tested using 



the overlap framework for smooth abelian subgroup configurations of SU(2), 
for lattice discretizations of instantons and for a cooled configuration. 

After the interest in the subject had been revived by the work of the 



overlap collaborators IT^, |l^, |4|, the last two years saw a wealth of 



of new articles. Low lying eigenvalues as well as complete spectra of the 



staggered fermion matrix were computed in |16, 17] and their distribution 
was successfully compared to the prediction of random matrix theory |18]. 
In |]19[ it was attempted to resolve the problem of exceptional configurations 
in quenched simulations by shifting the poles of the propagator due to zero 
modes. Studies of the behavior of fermionic observables in topologically 



non-trivial gauge field configurations can be found in [20, 21, 22 1. Results 
for the spectrum in classical configurations were reported in |2^, |2^ and 
in 24] studies in QED2 with dynamical fermions were performed. A 
detailed analysis of the lattice zero modes and their localization properties 
for the modified operators j^D is given in ]p^ ]. 

In this article we directly investigate properties of the unmodified Wilson- 
Dirac operator. In particular we compute complete spectra of the origi- 
nal Wilson-Dirac operator in sets of background configurations with known 
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topological charge. The features of the spectrum which are due to the topol- 
ogy of the background field are analyzed. We show the stability of these 
features under adding fluctuations to the gauge fields thus demonstrating 
their topological nature. The behavior of the spectrum in a fully quantized 
theory is discussed using QED2 as an example. 

The paper is organized as follows: In Section 2 we collect some known 
symmetry properties of the Wilson-Dirac operator and review the special 
role of the real eigenvalues. This is followed by Section 3 where we con- 
struct simple gauge field configurations to be used for numerically testing 
the lattice index theorem (Section 3.1) and discuss the numerics (Section 
3.2). In Section 4.1 we present our numerical results for smooth gauge field 
configurations, Section 4.2 contains the results for rough configurations and 
in Section 4.3 we perform a stability analysis of the topological features of 
the spectrum. Section 5 discusses the role of the index theorem in a fully 
quantized lattice gauge theory using QED2 as an example. The article closes 
with a discussion of the results and further perspectives (Section 6). 

2 Interpreting the spectrum of the Wilson-Dirac 
operator 

We work on a D-dimensional {D = 2, 4) lattice A with volume L^^^ x T. 
Lattice sites are denoted as x = (xi, ...xd)- The lattice spacing is set equal 
to 1. We write the fermion matrix M as M = 1 — kQ with hopping 
parameter k = (2m -|- 2D)^^, where m is the bare quark mass. Due to this 
simple relation between M and Q, in the following we will concentrate on 
the analysis of the hopping matrix Q. k is the hopping parameter, and the 
hopping matrix Q is defined as (spinor and color indices suppressed) 

D 

Q{x,y) = ^Qfi{x,y), 

For the 7^-matrices we use the chiral representation (see e.g. p^ ). The ma- 
trix anti-commuting with 7^, /x = 1, ...D is denoted as 75 for both D = 2, 4. 
In order to ensure reflection positivity, the fermions obey mixed boundary 
conditions, i.e. periodic in i-direction, i = 1, ...D — 1, and anti-periodic in 
D-direction. This is taken into account by the mixed periodic Kronecker 
delta lS.x±fi^y which is the usual D-dimensional Kronecker delta, but with 
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an extra minus sign for the links from xd = L to xd = ^- The gauge fields 
U^{x) are group elements (U(l) for D = 2 and SU(A^) for = 4) assigned to 
the links between nearest neighbors x,x + ft. They obey periodic boundary 
conditions. 

It is well known, that the hopping matrix Q is similar to its hermitian 
adjoint (note that Q is neither hermitian nor anti-hermitian) , and for 
even L and T also to —Q 

rsQFs = , EQE = -Q . (3) 

The matrices and H implementing the similarity transformations (^) are 
given by r5(x,y) = Sx^yj^, = {—l)^^~^'"^'^6x,y Both these matrices 

are unitary as well as hermitian and commute with each other. We stress 
that for odd L or T the second similarity transformation in (^) does not 
hold since the crucial condition Ax±fi,y{—l)^'^~^'"^^'^^'^^'"y^ = —Ax±fi.,y is 
violated for = 1 (or L,T). In order to be able to use the full symmetry 
(§) we work with even L and T in the following. The overall size of the 
matrix Q is r xr where r = L^^^ x T x D x N , with = 1 for gauge group 
U(l) (in = 2 only). 

Beside the symmetry transformations (^), Q also obeys a bound for its 
norm. With (|2|), using the fact that 1 + 7^ and 1 — 7^ are proportional to 
orthogonal projectors and taking into account the unitarity of the U^{x) it 
is straightforward to show |27], || Q/^ l|oo= sup||g||=i || Q^j^g ||= 2. The norm 



is defined to be the norm obtained by summing over all lattice, 
spinor and color indices and g is some (spinor) test function on the lattice. 
Using the triangle inequality one ends up with || Q ||oo < 21). This bound 
implies for an eigenvalue a of the hopping matrix Q (the corresponding 
normalized eigenvector is denoted as v) \a\ =|| Qv \\<\\ Q ||oo= 21). Thus 
the eigenvalues of Q are confined within a circle of radius 2D around the 
origin in the complex plane. 

It follows from the symmetries (|3|) that the spectrum of Q is symmetric 
with respect to reflection at both, the real and the imaginary axis. The 
eigenvalues of Q (and due to M = 1 — kQ also the eigenvalues of M) come 
in complex quadruplets or in real pairs. 



Let us now review the special role of the real eigenvalues and their eigen- 
vectors: Using the implementation of hermitian conjugation of Q as a sim- 
ilarity transformation one can show [|^] the following properties of the 
pseudoscalar matrix elements of eigenvectors Va,vp of Q with eigenvalues 
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a, f3 

Va^ ^5 Va 7^ Only for a G R , 

Va^ Vj3 7^ only for a = /3 . (4) 

only for real eigenvalues a. 

These equations should be compared to the chiral properties of the eigen- 
states of the (anti-hermitian) Dirac operator ^ in the continuum. Denoting 
an eigenstate with eigenvalue iE as ipE one finds (since 75 anti-commutes 
with ^) 

{ipE , 75 / only for E = , 

ii'E' , 75 i'E) + only for E' = -E . (5) 

When comparing @ and (|5|) it is immediately clear, that only the eigenvec- 
tors of Q (and thus of M) with real eigenvalues can be the lattice equivalents 
of the zero modes in the continuum. Based on this observation one can in- 
terpret the Atiyah-Singer index theorem on the lattice in the form [ p!o| ] 

v[U] = - R+ . (6) 

Here i?+ {R-) denotes the number of real eigenvalues a of the hopping 
matrix Q with positive (negative) matrix element vj^ Va in the physical 
branch of the spectrum, i.e. in the vicinity of +8 in the complex plane. 

Also the real eigenvalues in the other branches of the spectrum, corre- 
sponding to corners of the Brillouin zone different from (0,0,0,0), obey a 
simple pattern governed by the topological charge of the gauge field config- 
uration: The right hand side of (^) simply has to be multiplied by a factor 
corresponding to the degeneracy of the corresponding corner of the Brillouin 
zone (e.g. 4 for the (tt, 0, 0, 0)-type corners, 6 for the (vr, tt, 0, 0)-type corners 
etc.) and an extra minus sign for each vr in the coordinate of the corner 
(compare Table 1, Fig. |2| and the discussion in Section 4). @ is expected to 
hold for configurations with plaquette elements sufficiently close to 1, such 
that exceptional configurations are avoided. 

We remark that when analyzing the manifestation of the index theorem 
on the lattice the existence of Eq. (^ and its interpretation (P) makes the 
Wilson formulation more convenient compared to the staggered form of the 
fermion matrix. The latter formulation gives rise to an anti-hermitian matrix 
with only purely imaginary eigenvalues and thus lacks a natural criterion to 
identify the lattice equivalents of the zero modes. 
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Ui{x) 


= exp 




= exp ( 


U2{X) 


~ { exp 




~ { exp 



3 Numerical analysis of topological features of the 
spectrum 

3.1 Background configurations with non-trivial topology 

In order to test the form (|6|) of the Atiyah-Singer index theorem on the 
lattice, one has to compare the number of the real eigenvalues to the topo- 
logical charge of the gauge field configuration. For this enterprise it is ne- 
cessary to use a definition of the topological charge independent from the 
Dirac operator, such as the geometric definition ||2^, |2^. However, although 
conceptionally most satisfying, the geometric definition is rather costly to 
evaluate for Z) = 4 with gauge group SU(A^) |3^. For a first analysis of (P) 
in Z? = 4 we use sets of topologically non-trivial SU(2) gauge field configu- 
rations with known topological charge . These are configurations that are 
a lattice transcription of continuum SU(2) gauge fields with constant field 
strength on a torus. They are given by 

- ir [2ttj-^{x2 - 1) + ^i] 

-IT [27r-^(x4 - I) +^3. 

1 for X2 ^ L 

(ir [2^f(xi-l)+v92]) forx2 = L ' 

1 for X4 7^ T 

ir [27r|;(x3 - 1) + (/?4]) for X4 = T ' 

T is one (arbitrary) of the SU(2) generators (Pauli matrices). We remark 
that the spectrum of the hopping matrix Q does not depend on the choice 
of r, as can be seen from a hopping expansion of its characteristic poly- 
nomial. s,t are integers. We generalized the configurations given in |^] 
by introducing constant phases (p^ which leave the topological charge in- 
variant, but allow for a wider range of simple configurations where we can 
check the index theorem (^). For the configurations (^) the plaquette el- 
ements Ufj_iy{x) = Uf^{x)^U,y{x + e^)'^U^{x + ey)Uy{x) are constant with 
Ui2{n) = exp ( — i27rrs/L^), U^4,{n) = exp ( — i2-KTt/LT) and all other 
plaquettes U^u = 1. The topological charge is given by 

i^[U] = 2st . (8) 

The configurations (0) have action 4L^T[2 — cos(27rs/-L^) — cos(27rt/LT)] 
which for large L = T approaches the value 87r^(s^ + t^) > 87r^2|st| = 87r^|i^|, 
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thus obeying in this limit the lower bound for the (continuum) gauge field 
action in the sector v. 

The fields (^) are maximally smooth (in terms of their plaquettes). How- 
ever small fluctuations around the fields (0) should leave the topological 
charge and the properties of Q which depend on it (number of real eigenval- 
ues, chiral properties of the eigenvectors) invariant. For large enough L, T 
and small enough fluctuations the invariance of the (geometric) topologi- 



cal charge is an exact result: Liischer [28| gives the (non-optimal) bound 
Tr[l — Ufj_y{x)\ < 0.03 for a SU(2) configuration to be non-exceptional (i.e. 
the configuration can be assigned a topological charge; in the work of Phillips 



and Stone |2S] a similar bound is given). In order to change the topologi- 
cal sector, the gauge fields have to go through an exceptional configuration. 
Thus if the fluctuations are small enough, the topological charge remains 
invariant. 

In order to analyze configurations fluctuating around (0) we define new 
link variables 

U^,{x)rough = U^{x)^^^ U^{x)old , (9) 

where the smooth link variables U^{x)oid are given by (|^, and the fluctu- 
ations Ufj_{x)(^^-^ are SU(2) elements in the vicinity of 1, defined as (fij, j = 
1, 2, 3 denote the Pauli matrices) 

3 

U,{x)^,) = l2rf\x) + iY.'^J^^i^^ (10) 
where r\i^\x), j = 1, 2, 3 are sets of small random numbers with 
-e < r{f)(x) < e , j = 1,2,3 and rj^Hx) 



l-j2{rl^\x)y . (11) 
\ j=i 



The size of the fluctuations is bounded by e and in the limit e ^ we have 
f^/x(^)(e) ^ 1. In the numerical analysis described below, we will use values 
e = 0.15 and 0.3 and apply the roughening procedure (|9|) several times {rir 
times). Certainly the gauge fleld action increases with roughening. We 
found, that when applying several runs of the roughening procedure with 
different values of e, the average values of the action plotted as a function 
of e^/n^ lie on a universal curve. The roughening procedure thus behaves 
like a directed random walk towards higher values of the action. 

In order to avoid confusion, we finally remark that the configurations (0), 
(0) do not dominate the path integral of lattice gauge theory (see e.g. [p^. 
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32, 33]). They are merely a convenient tool to analyze the spectrum in gauge 



field configurations with known topological charge. 
3.2 Diagonalization of Q and checks 

Computing all eigenvalues of the fermion matrix and some of the eigenvec- 
tors is a challenging problem already for rather small, four-dimensional lat- 
tices. The fermion matrix for SU(A^) has size r xr with r = L^~^ xTxDxN . 
Studying e.g. SU(2) on a 4^ lattice means diagonalizing and computing 
eigenvectors for a 2048 x 2048 complex matrix. Storing the whole matrix in 
double precision requires already 64 MB of memory. Doubling the lattice 
size increases r by a factor of 16, and the required memory by a factor of 
256. Thus the memory requirements restrict one to rather small lattices. In 
particular we studied the spectrum for SU(2) gauge fields on lattices of size 
4^ and 4^ x 6. 

We remark that for hermitian matrices there exists a method with mod- 
est storage requirements: the Lanczos method without reorthogonalization 
p4[. It was successfully used by e.g. Kalkreuter and Simma |17] to di- 



agonalize the hermitian operator 75M. Unfortunately, for the unchanged, 
non-hermitian Wilson-Dirac matrix M (or Q respectively) the situation is 
less suitable for the Lanczos method without reorthogonalization. Barbour 
et al. ||9|, ^ report that they were not able to compute all the eigenvalues 
and they were forced to reorthogonalize. In that case, the Lanczos vectors 
have to be stored and storage requirements become comparable with stan- 
dard methods. Another typical problem for the Lanczos method without 
reorthogonalization is that the eigenvalues come with wrong multiplicities 
p4|, 17]. Since the exact number of degenerate eigenvalues is essential for 



our investigation, the use of standard routines for general complex matri- 
ces from the LAPACK package turned out to be the best choice: Q was 
first diagonalized (using LAPACK routines ZGEHRD, ZHSEQR) and then 
eigenvectors for particular eigenvalues were computed using ZHSEIN and 
ZUNMHR. Computing the spectrum (without eigenvectors) for the 4^ lat- 
tice takes typically 4 hours of computer-time on a MIPS RIOOOO CPU at 
180 MHz. The time for the computation of a single eigenvector to a given 
eigenvalue is only several seconds. 

Our investigation requires the distinction between real eigenvalues and 
eigenvalues with non- vanishing imaginary part. For some of the gauge field 
configurations we encountered complex eigenvalues with rather small imagi- 
nary parts of order 10~^ but there was never a problem distinguishing them 
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from truly real eigenvalues, which in the numerical computation have imagi- 
nary parts smaller than 10~^^ and do not have a complex conjugate partner. 
A second criterion to distinguish the eigenvalues is of course the pseudoscalar 
matrix element of the corresponding eigenvectors, which is of order 1 for real 
eigenvalues, but of order 10~^^ for eigenvalues with non- vanishing imaginary 
part. 

There are several ways to test the correctness and accuracy of the pro- 
grams. A first check we performed is comparing the numerical result for the 
trivial case (all link variables = 1) to the analytic result for the spectrum 
obtained by Fourier transformation. Our diagonalization routine reproduces 
the analytic result with an accuracy of 10^^^. 

A more advanced test makes use of results known from the hopping 
expansion (see e.g. H^). Using the fact, that 1 ± 7^ are proportional to 
projectors, one can show, that TrQ" is proportional to the sum over traces 
over products of links ordered along non-backtracking, closed paths of length 
n. In particular 

r r 

TrQ2 = ^CKi^ = , Tr = = ^ Re Tr ^7^ . (12) 

i=l 1=1 P 

When some of the sides of the lattice L or T have length 4, TrQ^ has some 
extra terms due to loops around the boundary. Sums over odd powers of 
eigenvalues Oi vanish trivially, since eigenvalues come always in pairs a and 
—a (compare (^)). We checked the quadratic and quartic sum rules, and 
found that they are obeyed with an accuracy of 10"^ for all spectra we 
analyzed. 

We remark, that from the hopping expansion it is also obvious that the 
spectrum of the Wilson-Dirac operator is gauge invariant: The characteristic 
polynomial det(Q — Q) has a hopping expansion in terms of TrQ". Each 
of these contributions is a sum over traces of closed loops of link variables 
and thus gauge invariant. This implies gauge invariance of the characteristic 
polynomial and thus of the spectrum. An alternative argument for the gauge 
invariance of the spectrum is to use the fact that the Wilson-Dirac operator 
transforms in the adjoint representation of the gauge group, and thus the 
eigenvalues remain invariant (the eigenvectors transform in the fundamental 
representation). 

4 Numerical results 
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4.1 Smooth configurations 

We start presenting our numerical results with the discussion of the spectrum 
for the smooth fields (^). We analyzed more than 50 of these configurations 
with different values of the parameter^ T,s,t,(pf^ for lattices with L = 4 
and T = 4, 6. Fig. |l| shows 4 spectra of Q on a 4^ lattice displaying different 
features of the spectrum which we will discuss below. 

Fig. ||.a shows the spectrum for the s = l,t = 1 configuration, which has 
topological charge = 2. The first thing to note are the expected symmetries 
of the spectrum with respect to reflection at real and imaginary axis. All 
eigenvalues lie within a circle of radius 8 around the origin in the complex 
plane, as discussed in Section 2 and their imaginary parts are bounded by 
4. The complex eigenvalues are well separated from the real axis, and many 
of them are degenerate (note that there are all together 2048 eigenvalues). 

There are 32 real eigenvalues, all of them multiply degenerate. They lie 
in the vicinity of the values -8, -4, 0, 4, 8. We denote these branches where 
the real eigenvalues are found as A*,B*,C,B,A as depicted in Fig. |2| and 
described in Table 1. As discussed in Section 2 Theorem, these branches 
correspond to the poles of the free propagator with periodic boundary condi- 
tions, situated at the corners of the Brillouin zone: (0, 0, 0, 0), (vr, 0, 0, 0) + 3 
permutations, etc. For the s = l,t = 1 configuration we find that the eigen- 
values in the branches A* , B* ,C, B , A are degenerate and the branches are 
occupied by 2, 8, 12, 8 and 2 real eigenvalues. The pseudoscalar matrix 
elements T^v have the same sign for all eigenvalues within a branch, with 
the signs distributed among the branches as — 1,-|-1,— 1,-|-1,— 1. All these 
features are as described in the lattice version (^) of the index theorem on 
the lattice and the right hand side terms of (^) for the principal branch are 
given by R+ = 0,R_ = 2. 

Fig. [l|.b and c show spectra for configurations with s = l,t = —3, having 
topological charge v = —6, and s = 1 and t = 0, with topological charge 
1/ = 0. For the former case we find that the numbers of real eigenvalues 
in each branch have increased to 6, 24, 36, 24, 6. The signs of the pseu- 
doscalar matrix elements are again equal for each branch, but now given by 
+1, -1, +1, -1, +1, since the topological charge is negative for this config- 
uration. Thus i?+ = 6, R- = 0. For the spectrum in Fig. |^.c, which is for 
a configuration with zero topological charge, we find no real eigenvalue and 
R-^- = R- = 0, compatible with (P). 

^ As remarked above, it can be shown that the spectrum of the Wilson-Dirac operator 
is independent of r, but different choices of r were used to check the program. 



10 



v = 2 
R^= 

• • • • " • •14 •••^ • • • 



-5 -4 -2 



v = -6 
R^= 6 

.•» 




-6 -4 -2 




v=0 v=2 
R= = 

' A - 



4 

1 




/? - 9 

1 


E 


-8 -6 -4 -2 

-1 

• t* 
• • • • • • 

l*\.;;.trtt: 

-4 


2 4 6 

!'•%- » •' 


8-8 ^ -6 -4^ -2 

-1 
-4 


2 ^4 6 ^8 

E 



Figure 1: Plots of the spectra of Q in the complex plane, for smooth 
background configurations (0) with s = l,t = 1 (a), s = l,t = —3 (b), 
s = l,t = (c) and for s = l,t = 1 (d). For Plots a, b, c all phases (p^ 
were chosen to be zero, while for Plot d we have 99^ = 27r/7, // = 1,2,3,4 
(small black dots) and compare it with the spectrum for the zero phase con- 
figuration (larger grey dots), u is the topological charge of the gauge field 
used (= 2st) and R+,R- are the numbers of real eigenvalues in the physical 
branch of the spectrum with positive and negative chirality. 
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Finally Fig. |^.d shows the spectrum for the s = 1, t = 1 case but now for 
all link variables we have chosen the non-zero value ip^ = 27r/7, /i = 1, 2, 3, 4 
for the phase in (^) (small black dots in Fig. |^.d). The large grey dots indi- 
cate the spectrum for the configuration with zero phase. The figure nicely 
demonstrates, that the phase can shift several of the complex eigenvalues, 
but the real part of the spectrum, which reflects the topological charge (un- 
changed by the phase) remains invariant. 

The general picture which emerges after our analysis of various spectra 
for smooth configurations (0) is summarized in Fig. ^ and Table 1 . It entirely 
confirms formula (^). 



A* B* C B A 




-8 -4 +4 +8 



Figure 2: Schematic picture of the branches containing the real part of the 
spectrum of Q. 



branch 


A* 


B* 


C 


B 


A 


R— — 


V 




Qv 




V 


corners 
of B.Z. 


(vr, TT, vr, vr) 


(vr,vr,vr,0) 
-|- 3 perm. 


(vr,^,0,0) 
-|- 5 perm. 


(^,0,0,0) 
-|- 3 perm. 


(0,0,0,0) 



Table 1: Properties of the branches containing the real part of the spectrum 
of Q. 



4.2 Rough configurations 

After testing the lattice version ^ of the index theorem for the smooth 
fields (0), we consider now the roughened fields (P). One expects, that for a 
small amount of roughening, which does not destroy the topological charge 
of the gauge field configuration, the relation (P) remains invariant. 

In Fig. ^ we show the spectra of Q (left-hand side plots) and the pseu- 
doscalar matrix elements v^V^v plotted at the positions of the corresponding 
real eigenvalues (right-hand side plots) for the s = l,t = 1 configuration 
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Figure 3: Behavior of the spectrum of Q (left-hand side plots) and the 
pseudoscalar matrix elements v'^T^v plotted at the positions of the corre- 
sponding real eigenvalues (right-hand side plots) for the s = l,t = 1 config- 
uration when the roughening procedure @) is applied. We show the results 
for roughening steps = 2,6 and 8, each with e = 0.3. For Ur- = 2,6 
(upper four plots) the total number of real eigenvalues in the spectrum is 
j^aii ^ j^aii _ 22 J according to (P), while for = S the gauge field becomes 
too rough, and (|6[) becomes violated giving rise to only 24 real eigenvalues 
in the spectrum. 
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having v = 2. We show the results for roughening steps = 2,6 and 
8, each step performed with e = 0.3 (see @). Comparing the spectra 
for the roughened configurations with the corresponding smooth spectrum 
Fig. m.a, one finds that the degeneracy of the eigenvalues is lifted. For 
rir = 2,6 the real eigenvalues are still more or less concentrated in the 
branches A* , B* ,C, B, A, but this structure becomes less pronounced with 
increasing number of roughening steps. For = 2,6 the total number of 
real eigenvalues is i?"'' + = 32 as expected from the index theorem 
(= 16 X For rij. = S the separation of the real spectrum into distinct 
branches has almost vanished, and the total number of real eigenvalues has 
dropped to i?"'' + i?^'' = 24. Obviously the configuration becomes too rough 
and (^) does no longer hold. It is also interesting to note, that the complex 
part of the spectrum approaches the real line. This trend continues when 
Ur is increased further, and the distribution of the eigenvalues becomes uni- 
form inside the disk of radius 4 around the origin for random configurations 
(compare |^). This observation nicely fits into the picture obtained from 
a strong coupling expansion. It can be shown, that the divergence of the 
pion propagator in strong coupling |36] is equivalent to a bound of 4 for the 
modulus of the eigenvalues of Q. 

The plots for the pseudoscalar matrix elements on the right hand side of 
Fig. ^ confirm the obtained picture. For = 2,6 the matrix elements v'^Qv 
all have the same sign within the branches of the real spectrum. For n,. = 8 
some of the real eigenvalues, and thus their matrix elements have vanished, 
and the remaining matrix elements do no longer obey the sign pattern of 
the smoother configurations. 

So far we presented only configurations where the matrix elements v'^T^v 
have the same sign for all real eigenvalues in a branch of the spectrum. We 
found, that our test configurations (|^, (^ also allow for cases where both 
signs occur within a branch. Due to the anti-periodic boundary conditions 
in time direction, configurations of the type s = 0, t 7^ 0, although they 
have = 0, give rise to spectra with an even number of real eigenvalues 
in each branch of the spectrum. Half of the corresponding eigenvectors 
have positive and half of them negative chirality and their contributions 
to the right hand side of (^) cancel. Fig. ^ shows the spectrum for the 
s = 0, t = —2 configuration after 1 roughening step (left-hand side plot) and 
the corresponding pseudoscalar matrix elements v'^T^v in all branches of the 
spectrum (right-hand side plot). Although the configuration has = there 
are 12 real eigenvalues in the spectrum. The index theorem on the lattice (^) 
is nevertheless fulfilled, since the signs of the pseudoscalar matrix elements 
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cancel in each branch of the spectrum. 
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Figure 4: Spectrum of Q (left-hand side plot) and the corresponding matrix 
elements v^T^v (right-hand side plot) for the s = 0, t = —2 configuration 
after one step of roughening with e = 0.3. The gauge field has u = and 
according to formula @, the sum of the signs of the matrix elements has 
to cancel for each branch of the spectrum. This feature is obvious from the 
right-hand side plot. 

The left-hand side plot in Figure Q also demonstrates nicely another 
important result of our analysis: The modulus of eigenvalues is no proper 
criterion to identify zero modes. From the spectrum plot in Fig. ^ it is 
obvious, that in the close vicinity of the real eigenvalues there cluster several 
complex eigenvalues. These complex eigenvalues have nothing to do with 
zero modes, since their pseudoscalar matrix elements vanish identically (~ 
10~^^ in our numerical analysis), while the matrix elements for the real 
eigenvalues are of order 1. 

We remark, that increasing the size of one of the lattice sides to 6 does 
not affect the structure of the spectrum. For all sufficiently smooth con- 
figurations we analyzed, the lattice version (^) of the index theorem was 
confirmed, clearly demonstrating that the structure of the real eigenvalues 
is the part of the spectrum which remains invariant within the topological 
sectors. 
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4.3 Stability analysis for the topological features of the spec- 
trum 

In this section we present results from the analysis of the spectra for a whole 
sample of roughened gauge fields. This serves to demonstrate that the topo- 
logical features of the spectrum are stable for a set of rough configurations 
(^) fluctuating around the smooth configurations 

In particular we analyzed the spectra for roughened gauge field configu- 
rations starting from the smooth configuration with s = t = l(z^ = 2). We 
then applied rir steps of the roughening procedure each with e = 0.15. For 
all together 15 values of in the range of = 16 to = 38 we computed 
spectra for 30 configurations at each n^. For each sample of 30 spectra 
(at fixed n^) we analyzed the topological features of the eigenvalues: The 
smooth configuration with t = s = 1 has all together 32 real eigenvalues, 
distributed among the branches as 2, 8, 12, 8, 2. The configuration has = 2, 
and thus according to (^) the quoted numbers are the minimal amount of 
real eigenvalues in each branch. In order to be compatible with (^) there is 
always the possibility that there are more real eigenvalues in a branch, but 
with different chirality of the corresponding eigenvectors, such that the dif- 
ference of left- and right-handed contributions gives the numbers 2,8,12,8,2 
for the branches. However, any spectrum with less than 32 real eigenvalues 
clearly violates (P). In almost all cases we analyzed, the breakdown of (P) 
was signalized by such a drop of the number of real eigenvalues below 32. 
For the following discussion we adopt as our criterion for a breakdown of 
(^) a violation of 

# of real eigenvalues = 32 . (13) 

We define the probability function p{nr) to be the probability of finding 
( p^ correct, and compute it for each sample of 30 spectra at fixed n^. 
Fig. ^ shows the probability p{ny.) (circles, scale on the left vertical axis) for 
different values of . In the figure we also display the average value of the 
gauge field action per plaquette S (diamonds, right vertical scale). 

The figure shows nicely, that there is a rather large number of roughening 
steps where the topological structure of the spectrum remains unperturbed, 
i.e. p = 1 (from = to = 22). During these roughening steps the 
action per plaquette S increases (starting with a value of 0.1 at n,. = 0) 
roughly by a factor of 25 which demonstrates the stability of the topolog- 
ical structure in more physical units. Then after 22 roughening steps the 
probability p for the index theorem (^ to hold drops sharply reaching a 
level of p ~ 0.05 at = 30. Under further roughening the probability 
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Figure 5: Probability p of finding (|^) correct (circles, left vertical scale) as a 
function of the number of roughening steps n^, each with e = 0.15. We also 
show the average value of the gauge field action per plaquette S (diamonds, 
right vertical scale). 



function keeps fluctuating around this value. The plot nicely demonstrates 
the stability of the topological structure of the spectrum. In particular real 
eigenvalues and the pseudoscalar matrix elements of their eigenvectors are 
seen to be rather stable which underlines their topological nature. 

We performed the same analysis (with less statistics) for roughened 
gauge fields starting from the smooth t = s = 2 configuration. Again 
we observed a large region where the topological structure of the spectrum 
remained stable, with a rather sharp drop of the probability function at the 
end of the region. We found that this drop was at a different value of the 
action per plaquette (approximately 25 % smaller), indicating that, at least 
for the lattice size we considered, there is no universal threshold in terms of 
the action where a topological sector gets destroyed. 
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5 Discussion of the overall picture in the fully quan- 
tized theory using QED2 as an example 

In the last section we have numerically investigated the form (|6|) of the 
index theorem on the lattice using test configurations for 4-D SU(2) lattice 
gauge theory. The set of gauge field configurations (0), we used is of 
course rather restricted. On the other hand, testing the obtained picture 
in a simulation of QCD4 is certainly too demanding for present computer 
technology. However it is possible to analyze the Atiyah-Singer index 
theorem (|6|) in the simpler model of QED2. 

QED2 is not only considerably cheaper to simulate, but also has a very 
simple expression for its topological charge: i^[U] = ^ J2x ^p(^)i where 
9p{x) = —ilnUi2{x) with 6p{x) G (— 7r,7r). Configurations where Ui2{x) = 
— 1 for one of the plaquettes are so-called exceptional configurations ||2^ , 
and no topological charge can be assigned to them. They are of measure 
zero in the path integral. 

Again it is possible to use the Atiyah-Singer index theorem in the con- 
tinuum as a guiding line of what to expect on the lattice. For QED2 there 
holds another index theorem, sometimes referred to as Vanishing Theorem 
[p8| . It states that there are either only left- or only right-handed zero 
modes. Furthermore, taking into account that in D = 2 there are 4 species 
of fermions, 1 in the physical branch and 3 doublers, one expects to find 

# of real eigenvalues = 4 . (14) 

We define p to be the probability of finding (U) correct and evaluate this 
probability for a whole ensemble of Monte Carlo configurations generated 



in a simulation of QED2 with dynamical fermions [24|. We computed p for 
several values of L = T and /?. The hopping parameter k was always chosen 
close to the critical value for the given L and /3. In Fig. |6| we plot^ the 
probability function p as a function of the dimensionless ratio (3 /L^ . 

It is obvious from this plot that when sending L ^ 00 and holding 
P/L"^ fixed (which corresponds to sending L ^ 00 while holding ^/L fixed), 
the probability function p approaches 1. This implies that (|T^ ) holds with 
probability 1 in the continuum limit. A detailed analysis shows, that 
this result is essentially uniform in v and also the chiral properties of the 
eigenvectors are conform with the expectations from Atiyah-Singer Index 
Theorem and Vanishing Theorem. 

^We thank A. van der Sijs for a discussion concerning the presentation of this data. 
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Figure 6: Probability p of finding (|T^ correct as a function of 13 /L?'. We 
show our results for lattice sizes L = 8, 12 and 16. The symbols are con- 
nected to guide the eye. 



It has to be stressed, that this result is not restricted to some subset of 
configurations (as is our analysis for QCD4) since it was established using the 
complete set of Monte-Carlo configurations. This demonstrates that the lat- 
tice version of the index theorem holds in fully quantized QED2. We remark, 
that in jl^] the interplay between the topological charge and the spectrum of 
the lattice Dirac operator was further investigated by demonstrating within 
the overlap formalism that the chiral condensate in the Schwinger model 
originates from the sectors with topological charge ±1. 



6 Summary and outlook 

The fact that for the Wilson-Dirac operator hermitian conjugation can be 
implemented as a similarity transformation with 75 has an important impli- 
cation for the chiral properties of the eigenvectors of the fermion matrix: 
Only eigenvectors with real eigenvalues have non-vanishing pseudoscalar 
matrix-elements and thus should be interpreted as the lattice equivalents 
of the continuum zero- modes. The lattice version of the Atiyah-Singer in- 
dex theorem thus connects the topological charge of the spectrum to the 
number of real eigenvalues. 
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In this contribution we analyzed those features of the spectrum of the 
Wilson-Dirac operator that are connected to non-trivial topology of the 
gauge field. In particular we demonstrated the stability of the real eigen- 
values and their corresponding pseudoscalar matrix elements under adding 
fluctuations to the gauge fields. This shows the topological nature of the real 
part of the spectrum. For the computationally less demanding case of U(l) 
gauge fields in two dimensions we analyzed the lattice index theorem for a 
complete set of configurations from a simulation with dynamical fermions. 
It was shown, that the [3 ^ oo limit is dominated by configurations where 
the lattice index theorem holds with probability one. 

The next logical step in a numerical analysis of the lattice index theo- 
rem would be to test the four-dimensional Wilson-Dirac operator with gauge 
field configurations from a Monte Carlo ensemble generated in a quenched 
or even unquenched simulation. On the analytic side it would be interesting 
to investigate the role of the topological part of the spectrum in fermionic 
correlation functions. In particular in a spectral decomposition of fermio- 
nic re-point functions one could try to separate the part corresponding to 
eigenvectors with real eigenvalues. The lattice index theorem (^) shows that 
these contributions are intimately connected to the topological sector of the 
background gauge field. Interesting insights into U(l) problem and Witten- 
Veneziano type formulas could be obtained. 

Acknowledgment: We thank Christian B. Lang for valuable discussions 
during the course of this work. 



Note added in the revised version: In the meantime five more con- 



tributions related to our work appeared as preprints. In [39| further studies 
of the spectral flow of the overlap hamiltonian are performed. In particular 
instanton configurations and gauge fields from unquenched simulations are 
considered. Furthermore the effect of an 0(a)-improving clover term is dis- 
cussed. [40 1 also analyzes the effect of the clover term on the level crossings 
for the overlap hamiltonian. Improvement seems to be a viable method to 
enhance the chiral properties of the eigenstates. In [^] the spectrum of the 
un-modified Wilson-Dirac operator is analyzed in the Schwinger model and 
perturbative arguments for the stability of the real eigenvalues are given. 
Finally in [42| the spectrum of the fixed point lattice Dirac operator in the 
Schwinger model is studied. 
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